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Abstract. Recent hexadecapole (vi) and elliptic (i^) flow measurements are 
used to constrain estimates for the degree of local equilibrium, mean free path A, 
and the viscosity to entropy density ratio (^) of the plasma produced in Au+Au 
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collisions at •y/sJvlv' = 200 GeV. The eccentricity-scaled flow coefficients 
and indicate that the plasma achieves a degree of local equilibrium within 
~ 5 — 10% of the value expected for a fluid with ^ equal to the conjectured 
lower bound of 1/47T. Estimates for A and 2 asa function of collision centrality 
and particle transverse momentum px, points to transverse expansion dynamics 
compatible with a strongly coupled low viscosity plasma. 

Keywords: elliptic flow, hexadecapole flow, viscosity, mean free path, local equi- 
librium 
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1. Introduction 

In central and mid-central Au+Au collisions at the Relativistic Heavy ion Col- 
lider (RHIC), matter is produced at energy densities well in excess of the value 
(~ 1 GeV/fm 3 ) required for a de-confinement transition to the quark gluon plasma 
(QGP). A strong indication that such a de-confinement transition indeed occurs, is 
the important role that quark-like degrees of freedom have been found to play in the 
transverse expansion dynamics leading to anisotropic flow. Such flow is routinely 
characterized by the even order Fourier coefficients 

= / e m( ^-* flp) \, n = 2,4 



where, 4> p is the azimuthal angle of an emitted particle, <&rp is the azimuth of the 
reaction plane and the brackets denote averaging over particles and events. The 
second (1)2) and fourth (v^) order coefficients characterize the magnitude of elliptic 
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and hexadecapolc flow respectively. 

At the highest RHIC collision energy of ^snn — 200 GeV, a universal scal- 
ing of elliptic flow, suggestive of constituent quark-like degrees of freedom in the 
collision zone, has been discovered [ [TJ [2] and is now well established for a broad 
range of collision centralities, particle species and transverse kinetic energies. For 
transverse momenta < 1.5 GeV/c, hydrodynamic calculations that model a lo- 
cally equilibrated QGP (with little or no viscosity) show good agreement with the 
data [ [3 01 [5] . A recent transport calculation which incorporates gluon dynamics [ 
|S] , also indicate good agreement with the magnitude and trend of centrality depen- 
dent i>2 data. Thus, the measurement and study of flow have been central to the 
confluence of experimental results which now bear evidence for the creation of the 
QGP in heavy ion collisions at RHIC. 

A less settled question which is still debated intensely, is the degree to which 
the QGP is thermalized [[7], and whether it is strongly or weakly coupled [HI [5]. 
Experimental constraints which allow estimates of the ratio of viscosity to entropy 
density 2 and the mean free path A, are crucial ingredients to the resolution of this 
question. In this contribution we show that the combined use of double differential 
t>2 and V4 measurements provide such constraints, and consequently, lend new and 
important insight to this question. 

2. Flow, local equilibrium & the coupling strength of the QGP 

The use of flow correlations as a probe for the nuclear equation of state (EOS) 
and the transport coefficients of hot and dense nuclear matter has been recognized 
for quite some time [ EJ QUI El 02] • The connection is made transparent in the 
framework of perfect fluid hydrodynamics where the conceptual link between the 
conservation laws (baryon number, and energy and momentum currents) and the 
fundamental properties of a fluid (its equation of state and transport coefficients) 
is straightforward. Much current effort to understand and determine transport 
properties are focused on several microscopic models (for recent reviews see for 
example Refs. [2J H3])- Hybrid approaches which involve the parametrization of 
deviations from hydrodynamic behavior (i.e full local equilibrium) are currently 
being studied as well [ [Pfl [23 [Mj ■ The latter exploits the fact that results from 
the Boltzmann equation reduce to those from perfect fluid hydrodynamics when the 
mean free path A becomes small [117]. 

2.1. Deviation from local equilibrium 

Further insight on the extent to which a system deviates from full local equilibrium 
can be obtained via simultaneous study of the scaling violations of V2 and U4. Here, 
the central idea is that partial equilibrium breaks the scale invariance of perfect 
fluid hydrodynamics (which requires full local equilibrium), and thus, gives rise 
to specific measurable scaling violations [Q3J [T5J [T5]. The quantification of such 
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scaling violations can then be used to constrain an estimate of the "degree" of local 
equilibrium and the transport coefficients. 

One such violation of the eccentricity-scaled second harmonic v 2 /£2, has been 
recently parametrized [[TH [15] via the Knudsen number K = \/R [119]. as; 

V2 = v% K-i 

£2 £2 K- 1 + Kq X ' 

where K~ x is proportional to the average number of collisions per particle N, in the 

. - . v h 

collision zone of mean transverse size R; A is the mean free path of these particles; 
is the eccentricity-scaled flow harmonic expected from perfect fluid hydrodynamics 
and Ko is a constant estimated to be 0.7 ± 0.03 with the aid of a transport model 
[ [16] . We have found that a modified form of Eq. [1] 



h r 



V2k = V 2k 
£2k £2k 



K~ 



i k 



k = l,2, ... , (2) 



provides a good estimate of the scaling violations for all even order harmonics 
investigated. 

Following the operational ansatz of Eq. [2] an estimate of the extent of the 
deviation from local equilibrium (degree of local equilibrium) can be obtained for a 
given centrality cut as: 

^ 1 fc = l,2,.., 



v\J (1+K/Ko) 

where K is obtained as a function of centrality via a fit to the eccentricity-scaled 
flow coefficients, as discussed below. 



2.2. Coupling strength and the ratio of viscosity to entropy density 

The shear viscosity r\ of the QGP medium reflects its ability to flow "freely" locally. 
This medium response to flow gradients is proportional to the range over which 
momentum can be readily transported transverse to the flow. Consequently, a 
necessary indication for a strongly coupled QGP would be the observation of a 
rather small value for the ratio of viscosity to entropy density 77/s. Here, it is 
noteworthy that this is a necessary but insufficient requirement because a weakly 
coupled system exhibiting anomalous viscosity could also give a relatively small 
77/s value [[S]. The ratio rj/s, of course, can not be arbitrarily small [HOI 121] 
because quantum mechanics limits the size of cross sections via unitarity. Therefore, 
operationally it is the extraction of a relatively small 77/s value in concert with a 
short mean free path (A), from data, that provides a robust indicator of a strongly 
coupled QGP. 

For a relativistic fluid, the ratio of viscosity to entropy density can be estimatd 

as: 

- w TAc s = KRTcs, (3) 
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where T and c s are the temperature and speed of sound respectively. Thus, the 
K values extracted from fits to eccentricity-scaled flow data, taken in concert with 
an estimate of T and a reliable EOS, can be used to evaluate ^. Similarly, the 
mean free path A = KR, can be estimated with the aid of the centrality dependent 
geometric value R. 

3. Proofing the extraction of Knudsen numbers from fits 
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Fig. 1. (a) Calculated values' of vs - N paTt for Several values of 2 

indicated. The dashed curves show the fits to these simulated data, (b) (- 
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0, 1, 2, 3))/(^(case — 1)) vs. A^ par t for for the fits performed in panel (a) see text. 
The dotted curves represent error bands. 

Reliable estimates for 2 ; \ and the degree of local equilibrium crucially depend 
on accurate extractions of the Knudsen number K , from fits to eccentricity-scaled 
flow data. Therefore, it is important to test the efficacy of the extraction procedure. 
Such a test is illustrated in Fig. Q] where we have performed fits to the eccentricity- 
scaled elliptic flow values ^ , recently calculated by Song and Heinz [ [T5] . The 
symbols in FiglDJa) show the |J values obtained from hydrodynamic calculations 
performed for 47r(^) values of 0, 1, 2 and 3 respectively (hereafter referred to as 
case — 0,1,2, and 3 respectively). The dashed lines show the corresponding fits 
obtained with Eq. [2] and the assumption [ [HI [22] that 
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Fig. 2. (a) «2 vs. -/Vpart and (b) V4 vs. -/V part for charged hadrons obtained with 
several reaction plane detectors for the pr selections indicated. The dashed lines in 
(c) and (d) show 10% and 20% error bands respectively. 



where a or (3 are fit parameters, S is the area of the collision zone, is the 
multiplicity density and N part is the number of participants. 

The curves in FigQJa) indicate a good fit to the simulated data for each case. 
We reiterate here that the two parameters of the fit are the scaled hydrodynamic 
limit ^ and and f3. The latter allows the determination of K for each centrality 

or value of iV par t. For each fit indicated in FigQTa), the value of ^ is within 5% of 
the calculated value of 0.23 (cf. solid circles in FigfT]). The dashed lines in FigOJb) 
show the ratio ( - (case — 0, 1 , 2, 3))/(- 2 (case — 1 ) ) determined with Eq. [3] and the K 
values obtained from the extracted values of 0. Within errors, the ratios shown in 
Fig. [ljb) are the same as the ratios of the -| values employed in the hydrodynamic 
calculations. Therefore, we interpret this agreement to be a good validation test of 
the reliability of the the extraction technique. 

4. Flow measurements and their scaling violations 

Reliable estimates for ^ , A and the degree of local equilibrium also demand robust 
measurements of the flow coefficients. Such measurements of v% and U4 for charged 
hadrons produced in Au+Au collisions at ^/snn = 200 GeV, have been recently 
carried out by the PHENIX collaboration. Fig. [2] shows a set of preliminary double 
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differential v-i and V4 data obtained from ~ 3.4 x 10 9 minimum-bias Au+Au events 
collected during the 2007 running period. These data were obtained via the reaction 
plane method of analysis [ [23] ; 

(cos(2A:(y>p - $ 2 ))) 
^ = (cos(2^ 2 -$KP))) fc = 1 ' 2 ' (4) 

where, ip p is the azimuthal angle of a charged track and $2 is the azimuth of 
the estimated second order reaction (event) plane. Five separate event planes 
were constructed with the aid of the PHENIX Beam-Beam Counters (BBC: 3.1 < 
I^bbcI < 3.9), Muon Piston Calorimeters (MPC: 3.1 < |?7mpcI <3.9), and the inner 

(i: 1.5 < |77RXNi| < 2.8), outer (o: 1.0 < |?7rxn | < 1.5) and combined (io) rings 
(North and South) of the newly installed PHENIX reaction plane detectors (RXN). 

The denominator of Eq.|3]is a resolution factor which corrects for the difference 
between the estimated $2 and the true azimuth $rp of the reaction plane [H3J. The 
three sub-event method [[24] was employed to obtain an estimate of this resolution 
factor (as a function of centrality) for each of the five event planes used in our 
analysis. For k = 1, the resolution factor for the combined reaction plane from 
both BBC's has an average of 0.33 over centrality, with a maximum « 0.42 in mid- 
central collisions [[53] [T]. The MPC and RXNi Q improve this resolution factor by 
about 35% and 100%, respectively. 

For mid-central collisions, a comparison of the double differential flow coeffi- 
cients V2,4,(pt, Ap ar t), shown for each event plane in Figs. Ufa) and (b), indicate 
excellent agreement (i.e much better than 5% and 10% for i> 2 and V4 respectively) 
over the broad range of pt selections shown. For very central and peripheral colli- 
sions, this agreement degrades to ~ 10% and 20% respectively. Here, it is important 
to note that; (i) the V4 signal is rather small (especially for low pt and central colli- 
sions) and is very difficult to measure precisely with a poor reaction plane resolution, 
(ii) The event plane resolution, for each event plane detector, achieves a maximum 
in mid-central collisions and worsens as one moves towards central and peripheral 
collisions. 

The rather good agreement between the measurements shown in Fig. [2] attests 
to their reliability and to the absence of a significant 77-dependent non-flow contri- 
bution. A prodigious non-flow contribution, such as from di-jets, would lead to a 
sizable difference in the V2 values obtained with event planes determined at different 
rapidity gaps (A77) with respect to the central arms [|2"5"]. 



4.1. Hydrodynamic scaling violations of V2 and 

To test for hydrodynamic scaling violations in the data, we divide the vi and V4 
measurements obtained with the RXNi event plan^l] (cf. Fig. [5]), by the respective 
eccentricity (ie. £2 and £4) for each centrality selection. Here, the guiding principle 
is the prediction from perfect fluid hydrodynamics that ^2,4 is proportional to the 



1 The RXNi event plane detector provides an optimal resolution across collision centralities. 
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Fig. 3. (a) v 2 /e2 vs. N palt and (b) 1)4/64 vs. N palt for several px selections as 
indicated. The dotted curves are fits to the data with EqJ5] 



initial spatial eccentricity £2,4 and is independent of the size of the collision zone R. 
Thus, perfect eccentricity scaling would be indicated by a flat dependence for both 
2* and ^ vs. AT part . 

For each centrality selection, the number of participant nucleons iV par t, was 
estimated via a Glauber-based model [US]. The corresponding transverse size R 
and £2,4 were estimated from the distribution of these nucleons in the transverse 
(x, y) plane via this same Monte-Carlo Glauber model [[26l [27], as well as via the 
factorized Kharzeev-Levin-Nardi (fKLN) [[28] model: 




(a 2 — a 2 ) 2 + 4<7 2 0^2 



,£ 4 =1- 



2 1 2 ' ^ 4 1 4 1 o 2 ' 

0- + a* 0% + cr* + 2a£ y 

where <7 X and <Tj, are the respective root-mean-square widths of the density distri- 
butions and a xy = xy — xy; here, bars denote a convolution with the density dis- 
tribution for a given configuration and averaging is performed over configurations. 
This procedure ensures that the fluctuation in the orientation of the initial almond- 
shaped collision zone is taken into account [[551 12Z]- For the Glauber calculations, 
the initial entropy profile in the transverse plane was assumed to be proportional 
to a linear combination of the number density of participants and binary collisions 
[[29]. This assures that the entropy density weighting is constrained by multiplicity 
measurements. 
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The eccentricity scaled V2 and i>4 values obtained with fKLN eccentricities, are 
shown in Fig. [3J For the lowest pr selection, they indicate relatively small scaling 
violations. However, these violations progressively increase as the px for charged 
hadrons is increased. That is, the data points slope progressively upward (from 
low to high iV par t) as the magnitude of the pr selection is increased. Relative to 
mid-central events, a large scaling violation of ^ is also apparent for the two most 
central bins ie. for N part > 250. We have traced this to a small overestimate of 
fluctuations which only impacts the value of £4 in the most central collisions as 
illustrated by the gray diamonds in Fig. [3J The latter shows the resulting values of 
^2 when an attempt is made to account for the overestimate by introducing a very 
small correlation between the £2 and £4 axes. 

A similar scaling test was performed with the Glauber model eccentricities. The 
resulting and |J values vs. N part , show trends which are relatively similar to the 
ones exhibited in Fig. [3J albeit with somewhat larger scaling violations. However, 
as discussed below, our fitting procedure provides a good constraint for choosing 
between the fKLN and Glauber model eccentricities. 



4-2. Estimation of the degree of local equilibrium, 2 and X 

As discussed earlier, the quantification of scaling violations is an important step 
for reliable estimates of -, A and the degree of local equilibrium. This has been 
accomplished by performing simultaneous fits to the eccentricity-scaled and V4 
data for each of the pr selections. For these fits, we follow the same fitting ansatz 
proofed in section [3j ie. we use Eq. [2j 



£2k £2k 



l\ 1 " ' 



K- 1 + Kg 1 



1,2 and K- 1 = /3N^ t . 



Here again, the scaled hydrodynamic limits and (3 are fit parameters. Note as 
well that j3 allows the extraction of K values as a function of i\r part , for each fit. 

The requirement of a simultaneous fit to the eccentricity-scaled data ensures 
that the same K values are extracted from the Vi and W4 measurements. Equally 
important is the fact that they provide a constraint for making a choice between the 
fKLN and Glauber eccentricities. That is, Glauber eccentricities result in relatively 
poor simultaneous fits while the fKLN eccentricities give very good simultaneous 
fits (-Radj ~ 1) as shown by the dotted curves in Fig.[3j These fits underscores the 
fact that data comparisons for different centralities, as well as ratios such as 
must take account of eccentricity differences. 

As pointed out earlier, each fit gives a value for the hydrodynamic limits — 

and -j-, and a (3 value which determines if as a function of N pSLIt . With these 
K values, the degree of local equilibrium and 2 can be estimated using Eqs. [H 
and [3] respectively. To estimate - at a given centrality, we assume a temperature 



Is the Quark Gluon plasma strongly coupled? 



9 




Fig. 4. (a) 



</E4 



0.5 1.0 1.5 2.0 
<p T > (GeV/c) 



100 



200 



300 



N 



part 



vs. pt, see text; (b) extracted bands for Air(rj/ s) vs. -/V part for 



severally selections as indicated; (c) extracted values of 4Tr(r]/s) vs. pt- 



T = 220 ± 20 MeV [ ;30J for the plasma when flow develops^, and use the the lattice 
EOS to estimate the associated value of c s — 0.47 ± 0.03 c. For any given centrality, 
A = RK is evaluated with the aid of the calculated transverse size R, obtained for 
that centrality. 

Figure U summarize the main results from our extractions with the fKLN eccen- 
tricities. It is noteworthy that the trends of the results obtained with the Glauber 
model eccentricities are essentially the same, albeit with different magnitudes. Panel 
(a) shows that, within errors, the ratio of the extracted hydro dynamic limits jqpjf/~p 

change little, if any, with hadron (pt)- The extracted bands for 4ir(r]/s) are shown 
as a function of N part in Fig.[4fb). The essentially flat dependence on AT part demon- 
strates the important role of the transverse size and the eccentricity of the collision 
zone. It further suggests that even though K increases as collisions become more 
peripheral, A does not have a strong dependence on collision centrality. 

In contrast to the iVp ar t dependence, Fig. [4jb) implies a relatively strong pr 
dependence for the extracted Airij]/ s) values, ie. i]/s increases by almost a factor 
of four over the indicated px range. This pt dependence is made more transparent 
in Fig. [3Jc) where we plot 47r(7y/s) vs. the mean value of the pt ranges indicated 
in Fig.QJb). The dashed-dot curve in the figure illustrates the quadratic nature of 
this pt dependence. Here, we wish to emphasize that this quadratic dependence on 

2 A centrality dependent T would results in an additional uncertainty on the extracted values. 
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Pt is not related to an intrinsic property of the QGP. Instead, it reflects the finding 
by Teaney [ [3T] that viscous corrections to ideal hydrodynamics grow as 

(£)'*■ 

Since our extracted K values include the pr-dependent factor given in the above 
expression, they give rise to the quadratic pr dependence of 4ir(r//s) observed in 
Figs-E^b) and (c). This also indicate that a pT-independent estimate of s) is 
obtained when the magnitude of pt ~ T. 

The slow rate of increase of 4ir(i]/s) at low pt (cf. Fig. BJc)) allows the use 
of the results for lowest pr hadrons to constrain the estimates 4tt(t)/s) = 1.3 ± 0.3 
and A = 0.25 — 0.3 fm for the plasma. The corresponding estimate for the degree of 
local equilibrium is within 5 - 10% of the value for a system having an 2 value equal 
to the conjectured lower bound of j-. To estimate the latter, we use the K values 
extracted from the fit to the simulated data (for 4ir(rj/s) = 1) shown in Fig. [TJ Our 
■q/s estimate is in good agreement with prior extractions [|4j [5], [6j HH [32J, [33]. Our 
A estimate also indicate that the low rj/s value is associated with a relatively short 
mean free path in the plasma. We interpret these observations as an important 
indication for a strongly coupled QGP. 

5. Conclusions 

To conclude, we have made detailed studies of possible hydrodynamic scaling viola- 
tions of the eccentricity scaled i>4 and V2 flow coefficients for charged hadrons. These 
studies validate the hydrodynamic scaling patterns expected for a nearly inviscid 
system close to thermal equilibrium ie. only 5-10% less than the value estimated 
for a system having an -q/s value equal to the conjectured lower bound of I/Att. 
They also provide the estimates 4ir(i]/s) = 1.3 ± 0.3 and A = 0.25 — 0.3 fm. These 
estimates suggest that the transverse expansion dynamics leading to anisotropic 
flow favor the creation of a strongly coupled QGP in collision zones for central and 
mid-central Au+Au collisions at ^/snn = 200 GeV. 
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